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Radiation from primary current sources above or at the interface of a homogeneous 
conducting dielectric half space has been studied extensively in connection with various 
communications problems. In many cases of practical interest, however, the surface bound- 
ing the conducting medium is irregular in configuration, the local deviation of the interface 
from the average being a function of position and time. The present paper describes a 
method for determining the effect of interface irregularities upon low-frequency electro- 
magnetic fields propagated within such a conductor. 

The electric field at a point within the medium is expressed as a stochastic integral over 
the bounding surface using the method of Stratton and Chu. As a consequence of certain 
restrictions which are principally geometric in character, a perturbation technique can be 
employed to simplify the integrand. Following the introduction of assumptions regarding 
t he st atisl ical properties of t he local surface deviations, approximate expressions are developed 
for the expectation of the subsurface field and the mean square deviation of the field ampli- 
tude. Numerical results for this latter quantity are obtained for the ease of a vertical 
electric dipole as the radiating source on the interface between free space and a moderately 
rough sea. It is shown that sea surface irregularities can produce a pronounced effect upon 
subsurface fields in some instances; a tentative physical interpretation of the Qumerical 
results is advanced. 



1. Introduction 

Propagation of electromagnetic waves above and 
within a homogeneous conducting dielectric half 
space due to certain prescribed elementary (time 
harmonic) current sources has been studied by a 
number of workers beginning with Sommerfeld [1909]. 
Modifications of SommerfelcTs original work by Weyl 
[1919] and Sommerfeld [1926], together with more 
recent extensions by Norton [1937], Bafios and 
Wesley [1954], Wait [1959a], Moore and Blair [1961], 
and many others, comprise a fairly exhaustive treat- 
ment of electric and magnetic dipole radiation in the 
presence of a semi-infinite conducting medium. 

In many situations of practical interest, the surface 
separating free space from the conducting medium is 
not smooth. In some instances it may be necessary 
to include local excursions of the surface from the 
mean in the derivation of expressions for field strengths 
both above and internal to the conductor. Even 
slight boundary irregularities will perturb the electro- 
magnetic field within the conductor from the value 
corresponding to a smooth surface. However, the 
comparative degree to which the fields are influenced 
by boundary roughness in various frequency ranges 
for a given conducting dielectric is a question which 
often remains to be answered. 
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A preliminary study of this problem has been 
reported by Wait [1959b] who considered a uniform 
plane wave in free space traveling in a direction 
parallel to the mean (plane) boundary of the con- 
ducting dielectric. The interface was taken to be a 
sinusoidal wave form with an amplitude much smaller 
than the period, with the magnitude of both being 
much less than the free space wavelength. Wait 
found that the resulting perturbation of the field 
within the conductor was proportional to the ampli- 
tude of the wavy surface in the first approximation. 

The present investigation is concerned with the 
propagation of electromagnetic waves within a 
conducting medium bounded by a rough interface 
in the presence of a prescribed current source. The 
precise configuration of the surface, however, is 
presumed to be unknown, with only its statistical 
properties being specified. As in the study 
performed by Wait, it is assumed that the vertical 
and horizontal scales of the surface irregularities 
are very small compared with the free space wave- 
length. A related problem of practical interest is 
the question of very low frequency signal propagation 
into the sea; for this reason, the conducting medium 
is referred to as the sea throughout the discussion. 
The general approach outlined here is applicable 
to other physical situations as well, provided certain 
relationships among the various parameters of 
importance are satisfied. 

For the low frequencies under consideration, sea 
water is characterized by an electrical conductivity 
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a s of about 4 mhos per meter and a relative per- 
mittivity e s of about 80. Conduction currents are 
much larger than displacement currents, and, as a 
result, the expression for the relative index of 
refraction N of the sea can be simplified to read 

where the time dependence of all field quantities is 
exp (—iwt). Cosequently, the effective propagation 
constant /5 can be written in the approximate form 



where 



p=k N=p exp(ify 



\a?€ / 



(2) 



(3) 



In terms of /3 , the effective wavelength X s in the sea 
and the skin depth 5 S can be expressed as 



and 



X s = 27rV / 2/3 - 1 = 1580/ A (/o m 
6 S = V2£ - 1 = 252/Vyom. 



(4) 
(5) 



Since the magnitudes of these parameters are of 
some importance to the developments which follow, 
plots are presented in figure 1 of X s and 8 S together 
with the free space wavelength X as functions of 
frequency /o over a wide frequency range. 




FREQUENCY (CYCLES PER SECOND) 

Figure 1. Variation of X s , 5 S and X with frequency. 

The geometry of the problem is illustrated in 
figure 2. A point source is located on the positive 
2 -axis and the observer is situated at the position 
(po, 0o 7 — kol). It will be assumed that the roughness 
of the ocean surface is unidimensional in character 
with the waves running in a direction which forms 
an angle <j> with the £ -axis. It will be convenient 
at a later stage in the analysis to employ a coordinate 



system whose origin is directly above the observer 
(at p , O , 0) with the ic-axis parallel to the wave 
direction. In this system the local height h of the 
ocean surface is expressible as a function of x alone. 
(The planes z = and 20=0 are coincident with the 
mean surface of the sea.) 

The electromagnetic field (or, alternatively, the 
Hertz potential) induced by a point source in the 
presence of a semi-infinite conducting dielectric with 
a smooth interface is usually given in a component 
form appropriate to a cylindrical coordinate system. 
In the present work, however, we shall be interested 
primarily in horizontal subsurface electric field com- 
ponents E xs and E ys ; i.e., the components in direc- 
tions parallel to and perpendicular to the wave 
direction, respectively. In the process of deriving 
expressions for these fields we shall have occasion to 
utilize the x, y, and z components of the free space 
electric field appropriate to the smooth surface prob- 
lem. Referring to figure 2, the transformation from 
cylindrical components (E p , Es , E Zq ) to Cartesian 
components (E XJ E y , E z ) is simply 



E x =E Pq cos 



-0)-^ o sin (00-tf) 



E y =E PQ sin (6 -<i>)+Ee Q cos (d -<t>) 
E Z =E Z „. 



(6) 



In section 2 approximate integral expressions are 
derived for the x and y components of the electric 
field within a moderately rough sea, in terms of the 
field (E x , E y , E z ) in free space induced by an arbitrary 
point source above a smooth sea. Subsequent to the 
introduction in section 3 of certain elementary 
assumptions regarding the statistics of the surface 
irregularities a quantity is derived which provides a 
measure of the effect of surface roughness upon the 
amplitude of the subsurface signal. Finally, in sec- 
tion 4, a specific source is assumed ; namely, a vertical 
electric dipole positioned at (x , y , z ) = (0, 0, 0+), 
transmitting at a frequency J in the VLF range. 



WAVE DIRECTION 




Figure 2. 



Geometry of the problem of propagation into^a 
rough sea. 
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Quantitative descriptions of the amplitude fluctua- 
tions of E Xi (the component of the subsurface electric 
field in the wave direction) are presented for varying 
degrees of surface roughness and various observer 
locations. 

2. Horizontal Electric Field Within the 
Rough Sea 

The electric field within the sea can be obtained 
by direct integration of Maxwell's equations using 
the method of Stratton and (Mm [1941], modified 
slightly to include the finite conductivity of the 
medium. It should be remarked that the derivation 
presented by these authors assumes a time depend- 
ence of the form exp (— icd) whereas, in our case, 
the changing structure of the sea surface produces 
a more complex temporal variation. For frequencies 
greater than a few tens of cycles per second, however, 
the intrinsic period of the electromagnetic "fields is 
far in excess of periods characterizing changes in the 
sea surface profile. Consequently, the integration of 
the field equations can be carried out with a harmonic 
time dependence as a "quasi-steady state' 7 approxi- 
mation which applies when the ocean surface profile 
is a slowly varying function of time. It will be 
understood in the sequel that the local deviation h 
of the surface from the mean plane 2=0 is a slowly 
varying function of time without explicitly exhibiting 
its dependence upon /. 

We shall employ the nondimensional variables: 

x'=x/\z \ 

y'=vl\zo\ 

z'=s/|so| 
h'=h/\z Q \ 



R'=R/\z \ = ^(x*+y*)/\z *\ 

r'=-z'x-?/y-(l+/Oz 



(7) 



where r' is the (normalized) vector from the point 
of integration to the observer. Moreover, tbe true 
vertical and horizontal scales of ocean surface 
roughness, <r and 8 , will be written as 



00=01201, 



-S\z \ 



(8) 



As a consequence of the nondimensionalization 
process, it will be convenient to employ the quantity 



B=B e 4=/3|s c 



(9) 



in place of p\z \ which appears throughout the 
analysis. 

Omitting the primary source terms, an expression 
can readily be derived for the electric field E s at the 
position of the observer in terms of the electromag- 
netic field vectors (E~, H") just below the ocean 



surface. Dropping the primes from the nondimen- 
sional variables defined in (7) the modified Stratton- 
Chu integral can be written in the form 



E, 



-vGX(nXE-) + (n-E-)vG]<IS. (10) 



In this equation, G is the solution of the three- 
dimensional Helmholtz equation appropriate to a 
medium characterized by a propagation constant 
B/\z \ and with a delta function source term : 



0- 



exp (iBr) 



It is easily shown that 



(lla) 



(lib) 



The ^-derivative of the local surface height h will be 
denoted by the symbol y. The local normal n to 
the surface can thus be written in the form (in the 
case of unidimensional roughness) 



n=(- 7 x+z)(l+y)- 1 /2. 



(12) 



The normal is directed outward from the medium, 
in the direction of free space. 

Certain simplifications can be effected by first 
utilizing the standard conditions of continuity of 
the tangential components of E~ and H" together 
with the conservation of charge flow across the 
boundary : 

nXE~ = nXE 



nXH=nXH 
n-E-=^n.E 



(13) 



where E and H are the electric and magnetic fields 
just above the surface S. Incorporating (lla), 
(lib), and (13) into (10), the expression for the 
electric field at the position of the observer appears 
in the form 



E 



= ^JX WoNo|(nXH) ~(^T)[ rX(nXE) 

-^(n-E)r]}^/S. (14) 

In order to extract quantitative information from 
(14) regarding the effect of a randomly rough ocean 
surface on the propagated signal, it is necessary to 
reduce the integrand to a combination of known 
deterministic functions and stochastic functional of 
the surface height h. To this end, we place the 
following restrictions on our considerations: 

(A) The sea surface profile is only moderately 
rough in the sense that the vertical scale a of the 
irregularity is at least an order of magnitude smaller 
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than the horizontal scale 5 (alternatively, the quan- 
tity y is to be considered at least an order of mag- 
nitude less than unity) ; at frequencies in the VLF 
range and lower, the scales of the surface irregulari- 
ties are much smaller than the free space wavelength ; 

(B) The horizontal displacement p of the observer 
from the primary source and the free space wave- 
length X are much greater than the skin depth 8 S ; 
consequently, the horizontal variation in amplitude 
and phase of the free space field 3 about £=0, y=0 
over distances a few skin depths in extent is negli- 
gible; 

(C) The location of the observer is such that the 
local excursions of the ocean surface from the mean 
can be considered small compared with p and \z \; 

(D) The depth of the observer is somewhat 
greater than a few 8 S so that 5 _1 is a small quantity 
compared with unity; and 

(E) The primary source frequency is sufficiently 
low to insure the smallness of l/|iV|. 

The upper limit of the frequency range is determined 
in part by restriction (E) and in part by the require- 
ment that the approximation expressed by (1) be 
valid. On the other hand, restriction (D) implies a 
lower limit on/ which depends upon the magnitude 
of |2 |. For example, for depths of the order of 50 m, 
the transmitter frequency should be greater than 
about 1 kc/s. When |2 | is considerably greater than 
this, the minimum frequency is correspondingly 
lower, extending to a few tens of cycles per second for 
depths in excess of 200 m or so. We also remark 
that the factor exp (iBr) in (14) provides for ex- 
ponential damping of contributions to the integral 
from portions of the integrand evaluated at x and y 
in excess of a few d s /\z \. This last observation, 
together with foregoing restrictions, implies that the 
quantities JS -1 , h, y, B 2 and iV -1 can all be regarded 
as at least first order in smallness. 

It is natural, therefore, to perform expansions 
within the integrand of (14) in series involving 
powers of the smallness parameters. For example, 
the inverse of the distance r has the expansion 



Correct to all orders, the normally directed differen- 
tial element of area ndS can be written as 



r- 1 =(l+2h+h 2 +R 2 )- 1/2 
= l-|i2 2 +|i2 4 - 



-A+l hB 2 +h 2 + 



(15) 



where terms of lower order than the second have 
been omitted. Simplification of the remaining in- 
verse powers of r proceeds in the same manner. 
The factor exp (iBr) may be expressed approxi- 
mately by retaining terms up to the first order in 
the expansion for r and then performing the substi- 
tution, giving 

exp (iBr) = exp \iB (l+i jR 2 +aYI. (16) 



3 In most cases of practical interest, the expression for the field at a height z 
involves quantities such as fcopo,/8po, and trigonometric functions of do. The in- 
equalities of restriction (B) insure the slow variation of the expression when po 
and po0o are altered by a few S$. An inspection of figure 1 reveals that the in- 
equalities are satisfied for po in excess of a few hundred meters over the frequency 
range of interest here. 



ndS= (—yx J r z)dxdy. 



(17) 



Moreover, we observe that according to restriction 
(B) the variation of the free space electric field in 
the horizontal direction is small throughout a region 
centered at (x,y) = (0,0) and a few 8 S in extent. 
This suggests that E(x, y, z=h) can be replaced ap- 
proximately by E(0, 0, z=h) which, in turn, can be 
expanded in the form 



E(0, 0, 2 =ft) ^*(A) = *(O)+i0'(O)+i A^"(0) + . . ., 

(18) 

where the prime indicates differentiation with respect 
to z. A similar approximation holds for the mag- 
netic field. Since the roughness scales of the sea 
surface are extremely small compared with the 
wavelength, the electromagnetic fields in free space 
are scarcely altered by the surface irregularities. 
Hence, the derivatives of <j>(z) can be calculated from 
the solution to the smooth sea problem with little 
error and the final result is a representation for E 
for substitution into (14) in which the only variable 
is h. The expression for the magnetic fielcTis treated 
in the same fashion. 

The approximate expression for E s is obtained by 
first substituting (15), (16), (17), and (18) into (14), 
together with the appropriate expansions of higher 
inverse powers of r and the series development of H 
which is the counterpart of (18). The final result 
for E s will be correct to the second order if those 
terms of the integrand are dropped which involve 
powers of smallness parameters which are lower than 
second order. 

The x- and ^-components of the subsurface elec- 
tric field are of principal interest here and can be 
written in the form (correct to the second order) of 
integrals oyer all x—y space in the manner just de- 
scribed. Since the random variable h depends only 
upon x, the integration over y can be performed at 
once. Some compactness of notation is achieved by 
replacing H with 



^=^H. 



(19) 



Omitting the algebraic details, the final result for 
the horizontal components of the subsurface electric 
field, following integration over y, can be expressed as 






11 

2 

fe=l 



*#* 2 



A A e*"e iBh Mx,h,y)dx, (20) 

%) — oo 

where j = x or y, and 
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and 



f 2 = x f Q =xy f^y^+iB-i) 

f^xt+iB- 1 f 7 =h 2 f u = x'+2iB- 1 x 2 -3B- 2 
fs=hy 

A xl =-xp y + (l + iB-')4> x 

A X2 =N-^ Z 

A x s=H,-(l+UB-^ x 

A x4 =^-^ v - (1 +2i^- 1 )0, + (l +iB- l W x 

A x5 =(l+iB-^ 2 

A x ,=0 

A x7 = — \p„ -\-\p' u — Wy +(t> x —4> x -f |<^' 

A xS =—(t>z+<f>'' z 

A xQ = -Wv+Wv+3<l>x-<t>x 

A x io= —<t> z 

A xn =<t>T—l^y 

Ay^h+il+iB- 1 )^ 

A y2 =0 
Ay^-h+x-il+UB-^ 

A y5 =f z 

Ay^-il+iB-^y 

A v1 =\l/ x —ip' x -\-Wx +4> u —4>y J rh ( i> , y 

Ay S =—\p Z +^' Z 

A V9 = §\I/ X — |^i+30 — <f)yy 

Ay 10= —\^ Z 

A v n=<i>y\-\tyx> 

In these expressions <f>] n) and ^j w) are understood 
to mean 4>j n) (0) and ^j tt) (0) where the differentiation 
precedes specialization of the argument. 

3. Sea Surface Structure and Electric Field 
Amplitude Fluctuation 

Field components E xs and E ys resulting from a 
particular wave structure could now be calculated 
simply by substituting the appropriate function for 
h into (20) and performing the remaining integrations 



numerically or in some other approximate way. 
In the present study the more general assumption 
is made that the sea profile is a stochastic function 
of x. A quantitative measure of the effect of bound- 
ary irregularity on the subsurface field can be 
obtained by calculating the mean square deviation 
(msd) of the field amplitude from its average value. 
Denoting mathematical (ensemble) averages (or 
expectations) by <g { } we have (after normalization) 

msd (\E Ja \) = g {IE*- g {E ia }] [E% 

-£*{Ei,}]}lg{E u }g*{E„} 
=[g [E js E%}/g {Ei a }g*{Ei,}]-l. (21) 

However, the analysis can proceed only after some 
assumptions are made concerning the statistics of 
the surface profile. Therefore, an elementary statis- 
tical model is introduced and certain averages 
necessary to the calculation of (21) are derived. 

The simplest classical model of ocean waves 
consists of an infinite succession of one-dimensional 
periodic progressive waves of nearly equal amplitude. 
In actuality, sea wave structure is almost never 
characterized by a single crest height or crest-to-crest 
distance. Instead, the profile of the surface through- 
out any local region is more realistically described 
as a superposition of a great many waves with a 
wide range of amplitudes and wavelengths. In the 
presence of moderate or high winds, the complexity 
of the wave configuration is enhanced and the 
surface loses nearly all appearance of regularity. 
The local departure of the sea surface from its 
mean value is therefore treated as a random variable 
in the present study. 

Sea surface structure is an extensive subject in 
itself and cannot be presented in any detail here 
(for an excellent summary the reader is referred to 
Defant [1961]). Suffice it to say that the statistics 
of ocean wave structure is dependent upon a great 
number of factors (including wind speed, duration, 
and fetch) and that the precise nature of this de- 
pendence is not thoroughly understood at the 
present time. For this reason (and for the sake of 
comparative mathematical simplicity) the surface 
height h will be described as a stationary Gaussian 
random process with a constant dispersion a 2 . 
The surface height has been assumed to vary with 
a single linear coordinate since swells and seas of 
moderate roughness often exhibit a general tendency 
toward unidimensional structure over fairly extensive 
regions. 

The restriction to "moderate" surface roughness, 

together with the assumption that 7=77 is a small 

quantity (for most realizations), requires that the 
variance (the expectation of the product h(x x )h(x 2 )) 
imply a value for the mean square local surface 
slope which is second order in smallness. It will be 
assumed here that the variance is a Gaussian function 
of the form 

g {hjio} =a 2 P (x h x 2 )=(T 2 exp [-fe-zOVS 2 ] (22) 
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where the subscript affixed to h signifies the point at 
which it is to be evaluated. It can be shown (see 
section 6, appendix A) that the mean square surface 
slope is given by 



g {y 2 } =lim <g {7172} =lim <? 



dxidxo 8 2 



(23) 



As a and 5 have been defined in this section, they are 
measures of boundary roughness in the vertical and 
horizontal directions and can be considered identical 
to the roughness scales previously employed. Hence, 
in accordance with the discussion in section 2, the 
mean square slope is a second order quantity with 
the variance defined as in (22). It should be noted 
that if an exponential variance were assumed, this 
condition would not be satisfied since the equalities 
in (23) would lead to an infinite value for the mean 
square slope. 

The calculation of the mean square deviation of 
electric field amplitude requires the expectations of 
both E js and the product E js E%. For the former 
quantity we have, from (20), 

A„ (" e'"'"S {«"*/.<*,*, -!■))<& (24) 

«y — 00 

and for the ensemble average of E js E* s , we find 
Bo 



<g {E js E%\ 



8?r 



- V2B0 



22 



+ ±jm/*-jn 1 1 ^ 


Bi 2 . B2 2 


iB x h x - 


- iBcJi 


2 


Jm\£\ 


>l,7l)/*fe>2 


T2)} 


(iX]( 


J'-j 


we have introduced the notation 








B X =B 










B 2 =B\ 









(25) 



(26) 



Since we are interested only in results correct to the 
second order, it is permissible to neglect lower order 
terms (such as m=8, n=4t) in the double sum in (25). 

The various expectations appearing in the inte- 
grands of (24) and (25) are readily calculated with 
the aid of the Karhunen representation (see e.g., 
Davenport and Root [1958]). The final results for 
those expectations required in the calculation of (24) 
and of (25) , subsequent to the deletion of terms which 
are lower than second order, are given in section 6, 
appendix A. It should be noted that the off-diagonal 
terms of the Cauchy product of E js and E% occur in 
complex conjugate pairs so that only one term of 
each pair need be calculated. 

With the aid of the equations listed in appendix 
A, the definitions following (20), and the relations 
given in (19), it is easily shown that the expectations 



in (24) have the simple form 



where 



g{e iBh j n (x,h,y)}=e **F m (z) (27) 

b=BW (28) 



and where F m {x) is a simple polynomial in x with 
complex coefficients. The integration over x can 
be performed at once, with the result 



&{E. 



JOS 



1 V2 R /V2 R b\ 



+ a 2 (l-B 2 a 2 )A j7 



+2~A J3 +iB^A H 



-2a 2 A. 



8 A 1- 



(29) 



The expectations appearing in (25) are somewhat 
more complicated; utilizing (19), (20), (22), and 
(26), as well as equations in section 6, appendix A, 
we find that 



g{e 



iB,h, 



' 4B * , Mx u h u ydftfaA, 72) } 

■^- Xi)W F m „(xi,x 2 ) (30) 



_ e &exp [ 



where F mn (x 1 ,x 2 ) is a polynomial in Xi and x 2 . Thus, 
the calculation of (25) has been reduced essentially 
to the problem of evaluating the integral 



00 

if 



1 2 Xl l 2 *%6exp[-Cz 2 -Zi) 2/ * 2 J 



F mn (x 1} X2)dx 1 dx 2 =Im n , say. (31) 



First making the change of variable 



——e v 



x 2 ~- 



u 



+ e 



V2 \ya 
and substituting into (31), we have 



(32) 



(32) 



S 2 

•\'2a 



//' 



where 



\'25 5 2 



F mn (u,v)dudv (33) 



(34) 



For realistic values of a (less than about 3 meters), 
the parameter b=Bla 2 is less than unity for the 
range of frequencies of interest here. Consequently, 
it is possible to expand the second exponential in 
(33) to obtain 
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=-*-**! ff 



■*■ mn /- 



e- u2 e- {a + l)c2 F mn (u,v)dudv. (35) 



Since 4 F mn (u, v) is a polynomial in u and v, it is a 
simple matter to calculate 4 /,„„. The final result for 
the expectation of E js E* is obtained by substituting 
the functions I mn into (25). This result, together 
with (29), provides all that is required for the calcu- 
lation of (21). The complete expression for the 
normalized mean square deviation of \E js \ is recorded 
in section 6, appendix B. In the next section a 
specific source is assumed and numerical results arc 
presented for various frequencies, sea states, and 
observer locations. 

4. Numerical Results for a Vertical Electric 
Dipole 

An application of the foregoing developments 
which may be of some practical interest is the case 
of a vertical electric dipole on the surface of the sea, 
transmitting at very low frequencies. It will be 
supposed that the observer's location corresponds 
to a horizontal displacement p which is several free 
space wavelengths in extent. The depth \z \ of the 
observer below the mean ocean surface is assumed 
sufficient to insure the smallness of /> _1 . 

The quantities of interest here are the mean 
square amplitude fluctuations of the horizontal com- 
ponents of the subsurface electric fields due to sea 
surface irregularities. Although at low radiofre- 
quencies the wavelength in free space is enormously 
Large compared with the scales a {) and 8 {) and the 
free space fields are hardly perturbed at all, we never- 
theless expect the surface structure to produce a 
pronounced effect upon the amplitudes of the sub- 
surface fields, particularly in the case of a vertical 
electric dipole. The reason for this can be seen by 
noting that, in the case of a perfectly smooth 
boundary, the dominant contribution to E xs and 
E ys arises from the horizontal component E p0 of 
the free space electric field, which is continuous 
across the surface z=Q. In the radiation zone, the 
2 -component of the electric field of a vertical 
electric dipole on the upper side of the boundary of 
a conducting dielectric is greater than E p0 by a factor 
of N, and this component is reduced by a factor of 
the order of 1/N 2 across the boundary of the conduc- 
tor. On the other hand, when the surface is charac- 
terized locally by a nonzero slope, the field E z0 will 
contribute to the local tangential electric field at 
the true surface. Since this component is coupled 
through the surface without change, it can contribute 
significantly to the horizontal component of the 
subsurface electric field. In the cast 4 of a vertical 
electric dipole over a unidimensionally rough sea, 
it is expected that this effect will manifest itself 
primarily in the ^-component of the subsurface 
field. This observation is verified by an inspection 
of (20), where the contributions to E xs from the 
projection of E z0 on the true surface are contained in 
the terms &=5, 8, and 10 which involve the product 
yE z . 



In order to perform numerical calculations, it was 
first necessary to obtain expansions of the type (18) 
for the electric and magnetic fields above a smooth 
conducting dielectric. The expressions used for (his 
purpose were those developed by Norton [1937] for 
a vertical electric dipole on a conducting half-plane. 
The coefficients of these expansions are required for 
calculation of the A jk recorded in (20). 

The results of several representative computations 
are shown in figures 3, 4, 5, and 6. These calcula- 
tions were performed to assess the effect of varying 
parameters on the mean square deviation of the 
amplitude of E xs . In figures 3, 4, and 5, the direc- 
tion of the ocean waves was assumed to be parallel 
to the radial displacement p of the observer from 
the dipole. The magnitude of p itself was fixed at 
500 km. Figures 3 and 4 reveal a steady decrease 
in the amplitude fluctuation with increasing depth 
and with increasing frequency for fixed sea states 
(constant o- n and 5 ). Figure 5 shows the variation 
of the fluctuation intensity with increasing correla- 
tion distance <5„ for several values of a . In contrast 
with figures 3 and 4 these curves are not monotonic, 
but instead indicate a maximum fluctuation level 
which depends upon the sea state, other parameters 
being fixed. 

Figure shows the effect of varying the angular 
coordinate O of the observer for fixed sea states when 
the ocean wave direction is parallel to the //,,-axis. 
The steady increase in the relative fluctuation level 
as the observer moves in the direction of decreasing 
() is not unexpected. When () is equal to zero, 
the observer is located on the ./' () -a\is. The ocean 
wave direction is <£ — 90° and therefore the ./--axis 
(and E xs ) is perpendicular to x . The source is a 
vertical electric dipole and for a smooth sea E xs 
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Figure 6. Variation of msd (|E XS |) with angular coordinate of 
observer. 

would be identically zero. In the present instance, 
however, the subsurface field E xs does not vanish 
because of the fact that the vertical component of 
the free space field can in general have a component 
tangential to the true surface. On the other hand, 
the mean field strength itself is at a minimum when 
# =0° and thus the relative fluctuation in intensity 
as defined here is greatest on the x -axis. Calcula- 
tions were also performed with p as the running 
variable, but only slight changes in relative fluctua- 
tion intensity were predicted. 

5. Physical Interpretation and 
Conclusions 

A few general observations can be made concerning 
the influence of ocean surface roughness on sub- 
surface field amplitude fluctuations by returning to 
figure 5. For a fixed rms surface wave height, the 
mean square deviation of \E XS \ is seen first to increase 
to a broad maximum and subsequently decrease as 
the horizontal scale of the surface roughness is in- 
creased. A possible explanation of this functional 
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behavior is provided by the considerations described 
below. 

Suppose the dipole to be transmitting at a fre- 
quency /o to an observer at a depth |2 |. On the 
average, the predominant contributions to the re- 
ceived signal arise from a roughly circular region of 
the mean sea surface. The boundary of this region 
is defined by the intersection of the mean surface 
with a cone whose vertex is at the observer's location 
and whose axis is colinear with the 2o-ax;is. The 
vertex angle a will be defined approximately by 
the relation 



COS a- 



N/l* 



(36) 



where \z x \ is the distance from the observer to those 
points on the mean surface where the contributions 
to the field at the observer have fallen off by the 
factor e~ l . In the case of a vertical electric dipole 
and nonnegligible surface slopes, the most important 
contributions originate in regions of maximum slope 
where the large vertical component of the free space 
field can contribute to E xs . Taking these contri- 
butions to be proportional to those of point sources 
in an infinite medium for order of magnitude esti- 
mates, the distance z x can be defined approximately 
as the solution of the transcendental equation 



_ t exp (— lgol/a a ) _ exp (— \zi\/8 s ) 



\H 



Zl 



(37) 



h=h sin (2wx/ch)). 



(39) 



where 8 S is the skin depth at frequency /<>. The 
radial extent a of the important contributing region 
of the sea surface is then 



«o=VNi| 2 — i^ol 2 - 



(38) 



Thus, at a fixed depth and frequency the radius a 
is determined. 

When 8 (the horizontal scale of irregularity, 
corresponding roughly to the length of a surface 
wave) is sufficiently small compared with a , the 
surface area of importance will contain a relatively 
large number of contributing regions, all radiating 
with more or less random phases. As 8 decreases, 
this state of affairs is enhanced and so also is the 
tendency toward destructive interference at the 
observer's location, resulting in a net decrease in 
the relative fluctuation in field intensity. On the 
other hand, if 8 is initially large compared with a 
and becomes increasingly larger, not only does the 
probability of finding regions of maximum slope 
within the area decrease but the excursions of the 
field at the observer due to these contributions are 
diminished since they are proportional to the deriv- 
ative of the local surface profile. Again, the net 
result is a decrease in relative fluctuation intensity. 

Between these two extremes is a transitional state 
in which, according to the foregoing remarks, the 
number of important radiating regions lying within 
the cone is of the order of unity. For the sake of 
discussion, assume the sea surface to be sinusoidally 
corrugated with a profile of the form 



Then a transitional state in which the fluctuation 
intensity is a maximum will be similar to the one 
illustrated in figure 7 where the observer is located 
below the region of maximum slope and where, for the 
sake of definiteness, one of the line segments \zi\ 
(a generator of the cone) is taken to pass through 
a point in the trough of the wave profile. The 
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Figure 7. Geometry of a representative transitional state. 

vertical and horizontal roughness scales <r and 
8 associated with the statistical model can be re- 
lated to the wave height h and wavelength d of 
the deterministic model by equating the ensemble 
averages of h 2 and y 2 with the mean values of these 
quantities as obtained from (39). Using (22) and 
(23) this procedure leads to the relations 4 



A -^V2 



0b 



do-^\2ir8 . 



(40) 



The validity of the foregoing observations can be 
tested in a simple manner. The transitional state 
which is illustrated in figure 7 should correspond 
roughly to a sea state for which the mean square 
deviation of lE^j is a maximum. The geometry in 
figure 7, together with (37), (38), and (40), is suffi- 
cient to determine the correlation distance (call it 
<5 trans which corresponds to this state when <r , 
|z |, and 8 S (or / ) are given. From figure 7 and 
(40) we have the relations 



do 



4a 



"V2tt V2tt|2 



(|so|— *d). 



(41) 



An approximate solution of (37) can be obtained by 
a single application of Newton's Method. Com- 
bining this result with (38) and (41), we obtain 



^n-MN-^J- 1 }''* < 42 > 



* Expressions such as these may be useful in relating statistical parameters to 
"apparent wave height" and "apparent wavelength" which oceanographers 
often associate with various sea states. 
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When plots such as those in figure 5 are drawn for 
several values of |2 | and/ > as weu as °"o> (with 8 
as the running variable) we should expect to find a 
maximum in each, located at a value of 8 which 
corresponds closely with 6 tran8 as given by (42). 
The dashed curve in figure 5 shows the locus of the 
maxima of curves corresponding to various a as 
calculated from (42). However, a test against a 
single variable is inconclusive. Calculations similar 
to those in figure 5 were performed for several values 
of / and 1 O I- Values of 5 trans were read from these 
curves and compared with calculations from (42). 
The results are presented in figures 8 and 9. The 
agreement in all three cases is seen to be satisfactory. 
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Figuke 9. Variation of transitional 5 with observer depth. 
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Figure 8. Variation of transitional 5 with frequency. 



In conclusion, it should be borne in mind that the 
quantity selected for evaluation in this paper, i.e., 
the mean square amplitude deviation of the horizon- 
tal components of the subsurface electric field, was 
chosen primarily to illustrate the method for assess- 
ing the effect of boundary roughness on subsurface 
signals. The same technique can be applied to any 



component of the electromagnetic field (or combina- 
tion of components) within the conducting medium. 
Moreover, the restriction to unidimensional rough- 
ness is not essential; bidimensional irregularities can 
be treated in precisely the same fashion. 
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6. Appendix A. Derivation of Certain Statis- 
tical Averages 

The surface height h(x) is assumed to be a separable 
random function, continuous for all x, with zero 
mean and variance o 2 p(x h x 2 ). In accordance with 
the Karhunen-Loeve Theorem, there exists an ex- 
pansion for h of the form 



in terms of the orthogonal process z p \ 

(r> 1 ZpZqJ ==0 PQ) 



(A-l) 



(A-2) 



where 8 pg is the Kronecker delta and where X" 1 and 
<f) p (x) are the eigenvalues and eigenf unctions of the 
integral equation 

J CO 
a 2 p(x h x 2 )<j> v {x 2 )dx 2 . (A-3) 
- 00 

An analytical procedure which is closely related to 
the one discussed by Hoffman [1955] provides us 
with the auxiliary relations (for a Gaussian Process) 



J iB 1 h 1 -iB 2 h 2 \ 



^Bl+Bl-ZB^p) _ 



= M(B U B 2 ) 



{ 



Zn€ 



>Bift, -iB 2 h 2 



■ =*X,[5i0 p (a!i)-5^ J ,(a; 2 ) JM(-Bi,J?a) 



—B 1 B 2 fj> p (x 1 )^> (I (x2) 
—BiBrf> p {x 2 )4 t (xi) 

+B 2 2 <t> p (x 2 )<l> q (x 2 )}}M(B 1 , B 2 ) 

(A-4) 

These equations, together with the definitions 

a 2 =2\lcj)l(x) (=const.) 

a 2 p(x h x 2 ) = 22\ p \ 9 (l) p (xi)<t) Q (x 2 ), (A-5) 

are sufficient for the computation of all expectations 
appearing in (25), subsequent to the deletion of 
terms of a lower order than the second. Averages 
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of the type appearing in (24) can be calculated after 
specializing (A-4) by setting B 2 equal to zero through- 
out. The various expectations required in the cal- 
culation of (24) and (25) are listed below. The 
results apply when h is a Gaussian process with zero 
mean and with dispersion and variance as given 
by(A-5): 

g{e im }=e^ BV =M{B) 

g {he im } =iBa 2 M(B) 

g{ T e iS "}=0 
g{h ye tB^ =0 

g {h 2 e iBh } =a 2 (l-5V)M(£) 

g { e «V>i-«%} = e -^ {Bl+Bl - 2B ' B2f, '=M(B 1 ,B 2 ) 
g \h,e iB ^- iB ^} =ic 2 (B 1 -B 2 p)M(B u B 2 ) 



<S {-ne*» 



-iB 2 <x 2 ^-M(B u B 2 ) 

4 <>P 



g {hme^-^h} =B 2 * i g {B.-B^MiB,, B 2 ) 
g {h 2 y 1 e tB ^- tI1 ^} 



dp 



(l+BtBtotp-BWWiBu B 2 ) 



S {h\e iB 



I"l- ! "2"2 



= <r 2 (l-Bl<7 2 + 2B 1 B 2 a 2 p-B 2 2 <r 2 p 2 )M(B 1 ,B 2 ) 

g {h 1 h 2 e iB ^- 1Ii2h *} 

= a 2 [p-B 2 yp+B 1 B 2 a 2 (l + p 2 )-ma 2 p}M(B u B 2 ) 



8 {7.72^ 



j«j — i£>2""2 > 



7. Appendix B. Normalized Mean Square 
Deviation of \E J0 \ 

The expression for the normalized mean square 
amplitude deviation of \E jX3 \ can be written in the 
form 

msd (\E j0 \) = {^ e~^ [g A jk A%P kk +2(Q m 

+ Qjil + $743 + Qj 52 + Qj&l + Qj 71 + Qm + Qj 111) 

-m+M>)\ m+M>)-\ (B-i) 
j 

The definitions of the A jk are recorded in (20). The 
appropriate expressions for the Q jmn are given by 



*&jmn J&eai \J±j m Jrlj n l mn ] 



(B-2) 



where the various P,„„ are defined as follows 

Fyi= Bl h ~a Ii+ 4 h+ M (a l2+h ) + B 2 /l 
P^Mi+BleVt+Byiu 



7^=2^7,-4^7,,- 



Blo 



g2 *" 



P«=~Y B A(l-i)Ii+(l+i)l7] 

P i3 =-<T 2 [(l-i)l2-Hi+(l+i)h 

+I 7 -(l+i)aI 3 +(l-i)aI 9 ] 

Ps2=-^B a% 
P tl =i^B^ 7 9 
P n =(< T *-iBl<r i )I 1 +2Blc i I 7 +iBlv*I 10 

P tl =i a i [(l+i)I 2 -a(l-i)I 3 +iI 1 

-{l-i)I s -a(l+ i)L-I 7 ] 

~(l-i)/ 8 +iJ|/,. (B-3) 
The I k appearing in these last equations are defined 

by 

Ii=I(a, b; 0, 0) 

I 2 =I(a, b; 1, 0) 

7 s =7(a, b; 0, 1) 

h=I{a, b;l,\) 

h=I(a, b; 2, 0) 

7 5 =7(a, b; 0, 2) 

I 7 =I(a+l,b;0,0) 

I a =I(a+l,b; 1,0) 

7 9 =7(a+l, 6;0, 1) 

7 10 =7(a+2, 6; 0, 0) 

7„=7(a+2, b; 0, 1) 
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where I(p, q; j, k) is the series 

/(p, q;j, k)=*D{j)D{k) ± 2- (p+r)-<*+*> (B-5) 



=o r! 



where 



z?o> 



(2j-l)l 
4*0-1)! 



D(0)=~ 



(B-6) 



Finally, the functions Lj and M, in (B-l) are ob- 
tained from 



2^ Bo 

Lj+iMj^e- 2 



[^, 1 +^- 3 (^) (l + i)-^(^ Boer*) (1-t) 

+4^(1-^0-^9(20+^11 (* J*)} ( B - 7 ) 
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